Unified description of hydrogen bonding and proton transfer by a two-state effective Hamiltonian 



Ross H. McKenzif] 

School of Mathematics and Physics, University of Queensland, Brisbane, 4072 Queensland, Australia 

(Dated: September 22, 201 1) 

An effective Hamiltonian is considered which describes hydrogen bonding and proton tranfer between two 
molecules due to the quantum mechanical interaction between the orbitals of the H-atom and of the donor (D) 
and acceptor (A) atoms in the molecules. The Hamiltonian acts on two diabatic states and has a simple chem- 
ically motivated form for its matrix elements. The model gives insight into the "H-bond puzzle", describes 
different classes of bonds, and empirical correlations between the donor-acceptor distance R and binding ener- 
gies, D-H bond lengths, and the softening of D-H vibrational frequencies. A key prediction is the UV photo- 
dissociation of H-bonded complexes via an excited electronic state with an exalted vibrational frequency. 



Hydrogen bonds play a key role in a diverse range of phe- 
nomena in physics, chemistry, molecular biology, and materi- 
als science. For example, hydrogen bonding is central to the 
unique properties of water, to protein folding and function, to 
proton transport fT|, to corrosion |2|, and to aspects of crys- 
tal engineering. Hydrogen bonds have properties that are dis- 
tinctly different from other chemical bonds. Many properties 
are poorly understood and recently the International Union of 
Pure and Applied Chemistry gave a new definition |4|, par- 
ticularly motivated by the occurrence of hydrogen bonding to 
a wide range of atoms other than oxygen (e.g., flourine, car- 
bon, nitrogen) and a wide range of bond lengths and energies. 
Gilli and Gilli have emphasized "the H-bond puzzle" which is 
stated as "the unique feature of the H-bond is that bonds made 
by the same donor-acceptor pair may display an extremely 
wide range of energies and geometries" |5|. A wide range of 
empirical correlations have been observed between different 
physical properties including bond lengths, binding energies, 
shifts in vibrational frequencies, vibrational absorption inten- 
sity, and NMR chemical shifts El 16). 

In this Letter, I consider a simple effective two state Hamil- 
tonian for a hydrogen bond between a donor denoted D-H and 
and an acceptor, which form a complex, denoted D-H- • • A. 
The Hamiltonian provides insight into the H-bond puzzle, 
gives a single description of different classes of H-bonds, and 
provides a semi-quantitative description of observed empir- 
ical correlations. The key independent variable (or physio- 
chemical "descriptor") is the distance between the donor and 
acceptor R (see Figure[T]i. 
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Figure 1. (Color online.) Definition of geometric variables for a 
hydrogen bond between a donor D and an acceptor A. 

Reduced Hilbert space for the effective Hamiltonian. Dia- 



batic states |8| (including valence bond states) have proven 
to be a powerful method of developing chemical concepts 
||9] [lOl . Previously it has been proposed that hydrogen bond- 
ing and hydrogen transfer reactions can be described by 
an Empirical Valence Bond model|ll-16| involving valence 
bond states. Here, I choose a basis consisting of two diabatic 
states which can be denoted as |D-H,A> and |D,H-A>. The 
latter represents a product state of the electronic states of a D 
atom and of a H-A bond in the absence of the D atom. These 
involve D-H and H-A bonds which have both covalent and 
ionic components, the relative weight of which depends on 
the distance r. The Morse potentials (see below) capture the 
associated energetics including the partial electrostatic char- 
acter of the H-bond. 

Effective Hamiltonian. The Hamiltonian for the diabatic 
states will have matrix elements which depend on the D-H 
bond length r, the donor-acceptor separation R, and the angle 
(f) which describes the deviation from linearity (compare Fig- 
ure [TJ. The functional dependence can be parametrised from 
quantum chemistry calculations lfT2l . However, the parame- 
ters may depend significantly on the level of theory used. The 
main point of this Letter is that one can obtain both a quali- 
tative and semi-quantitative description of hydrogen bonding 
using a simple and physically transparent parametrisation of 
these matrix elements. This approach highlights the quantum 
mechanical (covalent) character of the H-bond |7| and unifies 
H-bonding involving different atoms and weak, medium, and 
strong (symmetrical) H-bonds. The latter are sometimes char- 
acterised as covalent 4 electron, 3 centre bonds |9|. The model 
Hamiltonian has the advantage that it is straightforward to ex- 
tend it to describe nuclear quantum effects (including going 
beyond the Bom-Oppenheimer approximation), multiple H- 
bonds, and collective effects. 

The Morse potential describes the energy of a single bond 
within one of the molecules in the absence of the second (and 
thus the diabatic states). The two cases j = D,A denote the 
donor D-H bond and acceptor A-H bond, respectively. The 
Morse potential is 

Vjir) = L>j[exp(-2aj(r-roj))-2exp(-aj(r-roj))] (1) 

where Dj is the binding energy, r^j is the equilibrium bond 
length, and aj is the decay constant. The harmonic vibrational 
frequency lo is given by jiijp = ^Djoj where ^ is the reduced 
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mass. [For O-H bonds, lu ~ 3600 cm~^, D ~ 120 kcal/mol, 
a ~ 2.2/ A, tq ~ 0.96A.] A simple harmonic potential is 
not sufficient because the O-H bond is highly anharmonic and 
we will be interested in regimes where there is considerable 
stretching of the bonds. 

I take the effective Hamiltonian describing the two interact- 
ing diabatic states to have the form 



H 
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ADA{R,(b) VAir*) 



(2) 



where the diabatic states are coupled via the matrix element 

. s (R~^ r cos 6) 
Ada{R, 4>) = Ao cos(0)^ ^ exp(-&i?) (3) 



where r* = ^ R"^ + r"^ — 2rR cos ip is the length of the A-H 
bond (see Figure [T}, and h defines the decay of the matrix ele- 
ment with increasing R. This functional dependence on R and 
(f) can be justified from that for orbital overlap integrals if TTl 
together with a valence bond theory description of 4 electron 
3 orbital systems ||9l. With regard to the angular dependence 
I have assumed that the a overlap dominates the tt overlap. 
Roughly A^i^ is the overlap of the hybrid (s and p) orbitals 
on the D and A atoms and there will be some variation in the 
parameters Aq and h with the chemical identity of the atoms 
D and A. For the rest of the paper I focus on the case of linear 
H bonds {cj) = 0) and (|3) can be written as 



A{R) = Ai exp(-6(i? - Ri)) 



(4) 



where Ri is a reference distance, i?i = 2ro + 1/a. 

Potential energy surfaces. In the adiabatic limit the energy 
eigenvalues are 



E±[r,R) = -{VD{r) + VA{R-r)) 

±\{{VD{r)-VA{R-r)f 



-AA{Rf 



(5) 
■(6) 



I now focus on the symmetric case D = A and return to the 
asymmetric case briefly at the end of the paper. 

Figure |2] shows the adiabatic energy eigenvalues (poten- 
tial energy curves) E-{r,R) and E^{r,R) as a function of 
r, for three different fixed R values. Note three qualita- 
tively different curves, corresponding to weak, moderate, and 
strong hydrogen bonds. For moderate bonds (corresponding 
to i? ~ 2.6A for an O-H- • • O system) the energy barrier be- 
comes comparable to the zero point energy of the diabatic 
states, huj/2 ~ 1800 cm"^ ~ 0.04D. Hence, for smaller 
R quantum nuclear effects become important. The lowest 
curve in Figure |2] shows a ground state potential energy sur- 
face which has an extremely flat bottom, as is observed in 
quantum chemistry calculations for HsO^ and H5OJ ||T8!|. 

The value of Ai = O.AD ~ 2 eV used here for O-H- • • O 
systems is estimated from comparisons of the model predic- 
tions with experimental H-bond energies (see Figure|5]below) 
and vibrational frequencies (Figures |6] and |7]) |19|. This is a 
relatively large interaction between the diabatic states, and is 



comparable to values estimated from quantum chemical cal- 
culations at the MP2 level 1 12 |. The competition between this 
large energy scale and the large energies of stretched bonds is 
key to understanding H-bond properties. 




Figure 2. (Color online.) Potential energy curves for the diabatic and 
adiabatic states of a symmetric hydrogen bonded system. The hori- 
zontal axis is proportional to the difference between the length of the 
D-H bond, r, and its isolated value, ro. The vertical energy scale is 
D, the binding energy of an isolated D-H bond. All adiabatic curves 
are for Ai = 0.4D and b = a. The diabatic curves (dashed lines) are 
Morse potentials centred at r = ro and r* — R — ro and correspond 
to isolated D-H and H-A bonds, respectively. For parameters relevant 
to a O-H- ■ ■ O system the three sets of curves correspond (from top to 
bottom) to oxygen atom separations of i? = 2.9, 2.6, and 2.3 A, re- 
spectively, characteristic of weak, moderate (low barrier), and strong 
hydrogen bonds |5|. An important prediction of this model is the 
existence of the excited state with energy E+ (r, R) (upper curve). 



Insight into the H-bonding puzzle. Figures [3] |4j [5] show that 
varying R between 2. 4 A and 3.0 A can lead to a wide range 
of D-H bond lengths and hydrogen bond energies. In an ac- 
tual H-bonded complex the equilibrium R value will be deter- 
mined by the total energy which includes contributions from 
the potential surface shown in Figure |3] and other interactions 
not included in this model. Possible interactions include elec- 
trostatic interactions between the donor and acceptor, steric 
effects, and van der Waals interactions. In large molecules 
with intramolecular hydrogen bonding R may actually be de- 
termined largely by the skeletal geometry in which the donor 
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and acceptor atom are imbedded. Similarly, in water-hydroxyl 
overlayers on metal surfaces, the distance R is largely deter- 
mined by the lattice constant of the substrate IJ). This pro- 
vides insight into the origin of the H-bond puzzle because 
the identity of the donor and acceptor atoms determines the 
Hamiltonian parameters except R, whose equilibrium value 
will be determined by residual interactions. 




Figure 3. (Color online.) Contour plot of the ground state potential 
energy surface (for a symmetric donor acceptor-system) as a function 
of the D-H bond length a{r — ro) (horizontal axis) and the donor- 
acceptor distance a{R — 2ro) (vertical axis). Note that as R varies 
a wide range of equilibrium bond lengths rm{R) are possible. For 
this plot Ai = OAD and b — a. The contour spacing is 0.07D and 
darker shades represent lower energies. 
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Figure 4. (Color online.) Correlation between the H-bond energy 
and donor-acceptor distance R (in A). Solid line is for the model 
Hamiltonian with b = a = 2.2/ A, Ai = OAD. The dashed curve is 
an empirical relation, Ehb = E^ax exp(fc(i?o — R)) with fc = 5.1, 
Ro = 2.4A, and Emax = 27 kca^mol = 0.25D for O-H- ■ ■ O 
bondskS). 

H-bond energies. Figure |4] shows a favourable comparison 
between the calculated binding energies as a function of R 
with an empirical relation for O-H- • • O bonds |5 1. 

Bond lengths. In the adiabatic limit, the equilibrium D-H 
bond length r for a fixed R, is determined by the minimum 
of E_ (r, R) (see Figure[5|l. With decreasing R the D-H bond 
length increases from its non-interacting value tq ~ 0.96A 
to R/2 for a symmetric bond associated with a barrierless 



potential well. Figure |5] shows that the bond lengths calcu- 
lated from the minimum of (r, R) are significantly less 
than those observed experimentally for R < 2. 6 A. This dif- 
ference is attributed to the importance of quantum nuclear ef- 
fects which become when the energy barrier is comparable to 
the zero point energy of the D-H stretch vibration. A signa- 
ture of such quantum effects are isotope effects. Indeed this 
can be seen by comparing the crystal structure of CrH02 and 
CrD02; in the former the O-H-O bond appears to be sym- 
metric and R = 2.49 ± 0.02A, whereas the O-D-O bond is 
asymmetric with an O-D bond length of 0.96 ± 0.04 A, with 
= 2.55 ± 0.021 lEOl. 
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Figure 5. (Color online.) Correlation between the D-H bond length 
r and the D-A distance R. Both lengths are in units of A. The solid 
curves is the bond length deduced from the minimum of the adiabatic 
potential for b — a and A(i?) = 0.4D. For the moderate to strong 
H-bonds which occur for R < 2. 5 A quantum nuclear motion will 
significantly increase the D-H bond length because the energy barrier 
becomes comparable to the zero point energy. The dots are experi- 
mental data for O-H- ■ ■ O bonds in a wide range of crystal structures 
and are taken from Figure 6 in 1221 . The dashed line corresponds to 
symmetric bonds (r = R/2). 

Vibrational frequencies. As the distance R decreases there 
is a significant softening of the frequency of the D-H stretch 
vibration. In the adiabatic limit this frequency is given by the 
curvature at the bottom of the potential E- (r, R) (see Figure 
[6|. The latter has been proposed as a measure of the strength 
of an H-bond (3). Generally, it is expected that when the quan- 
tum nuclear motion is taken into account the actual vibrational 
frequency will be less than adiabatic harmonic frequency. 

A key prediction of the model. The energy eigenvalue 
E+ (r, R) describes the potential energy curve of a low ly- 
ing excited state (see Figure |2]) which reflects the quantum 
mechanical character of the H-bond. This state would not ex- 
ist if the H-bond is purely classical and electrostatic. This 
"twin excited state" is the analogue of the ^B2u state in ben- 
zene, of the state in semibulvalene |9|, and of the low- 
lying electronic state in the Creutz-Taube ion associated with 
delocalized mixed valence 1211 . Hence, for strong hydrogen 
bond complexes there should be an electronic excited state 
with energy of approximately 2A ^ 4 eV (corresponding to a 
wavelength of about 300 nm). The transition dipole moment 
is approximately equal to twice the ground state dipole mo- 
ment of an isolated OH bond, 2doH ^ 2 Debye, suggested a 
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Figure 6. (Color online.) Softening of the D-H stretch frequency Q 
(in cm~^) with decreasing donor-acceptor distance R (in A). Solid 
line is the harmonic frequency for the model Hamiltonian with b = a 
and Ai = OAD. The dots are experimental data for a wide range of 
complexes and are taken from Figure 4 in Ref . (6] 



Figure 7. (Color online.) Softening of the D-H stretch frequency 
Q (in units of cm~^) with decreasing difference between the proton 
affinity of the acceptor and the donor (e in units of kJ/mol). The 
solid line is the harmonic frequency of the model Hamiltonian with 
R = 2ro + 1/a ~ 2.4A, Ai = OAD a. Da = 120 kcal/mol. The 
dots are experimental data 1231 . 



significant absorption intensity. The curve {r, R) has only 
a single minimum as a function of r and so in this excited 
state the H-atom will be delocalised between the donor and 
acceptor For strong bonds the curvature of the adiabatic po- 
tential of the excited state is clearly larger than that of the 
ground state. Hence, the corresponding vibrational frequency 
will be larger than in the ground state, as it is in benzene and 
semibulvalene [gj. The intensity of the transition will be re- 
duced by Franck-Condon factors describing the overlap of the 
nuclear wave functions in the ground and excited state. For 
weak bonds these overlaps may be small. Since A(i?) de- 
creases with increasing R the energy of this excited state will 
also decrease with increasing R, implying that exciting to this 
state will lead to photo-dissociation of the H-bond. 1 am un- 
aware of any experimental investigations that have looked for 
this excited state in the UV. 

Asymmetric bonds. We now consider the case where the 
proton affinity (PA) of the donor (Do) and of the acceptor 
(Da) are unequal, i.e., e = Do - Da 7^ 0- GilH and Gilli (5] 
have noted the "PA, pKa equalisation principle": the strongest 
hydrogen bonds occur when the donor and acceptor have the 
same proton affinity. This occurs naturally in our two-state 
Hamiltonian. For example, for very short D-A distances and 
large A, the H-bond energy is approximately, e — \/e^ + A^, 
which has its largest value for e = 0. Figure [7] shows the soft- 
ening of the D-H stretch frequency as a function of e, com- 
pared to experimental data ll23l . 

Quantum and isotope effects. The above discussion treated 
the nuclear degrees of freedom classically, but noted that 
quantum nuclear effects may have a significant effect on equi- 
librium bond lengths, particular for strong H-bonds. Signifi- 
cant isotope effects are observed expeiimentally [6] and arise 
due to the zero point motion of the hydrogen atom. Further- 
more, quantum nuclear effects may play an important role in 
water L24il . in water-hydroxyl overlayers on metal surfaces 
(|2l, and in some proton transfer reactions in enzymes ||25l . 
For strong bonds where the potential is very anharmonic quan- 
tum nuclear effects will also be important |3|. The model 



Hamiltonian provides a natural means to describe these effects 
if the the hydrogen atom co-ordinate r is treated quantum me- 
chanically. The harmonic limit corresponds to a spin-boson 
model which has an analytical solution in terms of continued 
fractions [27 1 . The fully quantum Morse potential has an exact 
analytical solution and an algebraic representation in terms of 
creation and annihilation operators p28 |. Hence, an algebraic 
treatment of the quantum version of the model Hamiltonian 
may also be possible, because the off-diagonal terms are inde- 
pendent of r. 

In conclusion, the relatively simple effective Hamilto- 
nian considered here provides a unified picture of a range 
of phenomena associated with hydrogen bonding. Futher- 
more, it predicts an excited state which should lead to 
photo-disassociation of an H-bonded complex. Future work 
could consider non-linear bonds and the associated vibrational 
bending modes, and correlations between R and vibrational 
absorption intensities and NMR chemical shifts. It would be 
worthwhile to provide a more rigorous justification of the dia- 
batic state Hamiltonian from quantum chemistry. This may be 
done in a similar systematic manner as has been done for the 
excited states of the chromophore of the Green Fluorescent 
Protein f26]. 
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